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Abstract 

An array of resistively and capacitively shunted Josephson junctions with nonsinusoidal current- 
phase relation is considered for modelling the transition in high-T c superconductors. The emer- 
gence of higher harmonics, besides the simple sinusoid I c sin cj>, is expected for dominant d-w&ve 
symmetry of the Cooper pairs, random distribution of potential drops, dirty grains, or nonstation- 
ary conditions. We show that additional cosine and sine terms act respectively by modulating the 
global resistance and by changing the Josephson coupling of the mixed superconductive-normal 
states. 

First, the approach is applied to simulate the transition in disordered granular superconductors 
with the weak-links characterized by nonsinusoidal current-phase relation. In granular supercon- 
ductors, the emergence of higher-order harmonics affects the slope of the transition. Then, arrays 
of intrinsic Josephson junctions, naturally formed by the Cu02 planes in cuprates, are considered. 
The critical temperature suppression, observed at values of hole doping close to p = 1/8, is inves- 
tigated. Such suppression, related to the sign change and modulation of the Josephson coupling 
across the array, is quantified in terms of the intensities of the first and second sinusoids of the 
current-phase relation. Applications are envisaged for the design and control of quantum devices 
based on stacks of intrinsic Josephson junctions. 
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I. INTRODUCTION 



Arrays of Josephson junctions are under intensive investigation for potential implemen- 
tation as quantum bits and for modelling several phenomena in superconductive films^— . 
In particular, the resistive transition has been described by means of resistively and capac- 
itively shunted Josephson junctions arrays (RSCJ)^~— . In the conventional RSCJ model, 
the Josephson current is the simple sinusoid Is{<f>) = -f c sin<^>, where I c is the critical current 
and <f) = 02 — Oi the phase difference of the superconductor order parameters Ai exp(z$i) and 
A 2 exp(z6 , 2 ) . Sign and magnitude of I c are affected by the gap function symmetry and rel- 
ative orientation of the superconductor electrodes. For conventional phase-coherent pairing 
with s-wave symmetry and A x = A 2 = A, the critical current is given by the Ambegaokar- 
Baratoff expression I c = ir A 2 /2e-R tanh[A/2£;T], with R Q the normal-state resistance^. For 
unconventional high-T c superconductors, the internal structure of Cooper pairs most likely 
agrees with predominant d-w&ve symmetry that might originate deviations in the supercon- 
ductive and normal branches of the current- voltage characteristics. Moreover, ferromagnetic 
impurities, grain boundaries and interfaces, vortex cores, impurities and far-from equilib- 
rium conditions cause the onset of higher harmonics^— . The current-phase relation is given 



with f(E) the electron energy distribution and Im[/g;(^)] the spectral current, which depend 
on material, geometry and nonequilibrium conditions. Eq. ([T]) can be written as an n-order 
Fourier series^— : 



When the sum is restricted to n — 1, I n sm(n<j)) reduces to the familiar sinusoidal Josephson 
current 7 c sin</>. J n cos(n<f)) is the quasi-particle-pair-interference current, vanishing when 
the pair-symmetry is not broken. For s-wave superconductors and n — 1, I n shows a loga- 
rithmic divergence at V = 2A/e, whereas J n is zero for V < 2A/e with a discontinuity at 
V = 2A/e both in normal (0) and ferromagnetic (it) junctions^^— . It has been estab- 
lished that higher harmonics are important in cuprates. At the same time, non-monotonic 
temperature dependence of Josephson current also appears in d-wave system^^— . For 



by: 
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unconventional superconductors with prevalent ti-wave pairing, the harmonic sin 20 is criti- 
cally enhanced by the presence of midgap Andreev resonant state^i^ and can even dominate 
over sin0, as found Deviations from the sinusoidal shape can be more easily 

observed at temperatures below T c because, in general, these effects are of the second order. 
Nonetheless, in the vicinity of T c , they have been observed in normal-metal weak-links, as 
a consequence of the depairing either by proximity effect by supercurrent or in long junc- 
tions or in nonequilibrium conditions^. Spin-singlet/spin-triplet superconductor—^ and 
superconductor/ferromagnet hybrid structures have triggered considerable interest in recent 
years for their potential spintronics applications as they allow for tuning the critical current 
via the electron spin. If the metal between conventional superconductors is magnetic, the 
symmetry is broken and the current takes the more general form: I s = 7 sin(</> + <ft ). The 
phase-shift 0o is proportional to the magnetic moment perpendicular to the potential of the 
spin-orbit coupling^— . 

In this work, a model of the superconductive-resistive transition based on a network of re- 
sistively and capacitively shunted nonsinusoidal Josephson junctions (RSCNJ) is considered. 
Such a network could be relevant when the overall and concurring effects of above described 
phenomena should be taken into account. The ultimate scope being the consistent descrip- 
tion of several experimental evidences, that cannot be accounted for by the simple sinusoidal 
coupling. The appropriateness of the resistively-capacitively picture in the presence of the 
nonsinusoidal current-phase relies on the occurrence of the macroscopic quantum tunnelling 
in high-T c materials with <i-wave symmetry, whose experimental evidence has been reported 
only very recently^ 1 ^. 

Arrays of weak-links in polycrystalline superconductors and intrinsic Josephson junctions 
in cuprates are considered as prominent examples. 

In polycrystalline superconductors in the vicinity of the transition, nonequilibrium ef- 
fects make the relevant properties of the weak-links spatially and temporally dependent 
on the external drive^^— . Hence, when a polycrystalline superconductor undergoes the 
transition, the onset of higher harmonics may occur according to the local voltage, geom- 
etry and chemistry of the grains. The pair-interference current J n cos(n0) emerges when 
the pair-symmetry is broken and comes into play when the junctions are partly dissipative. 
Therefore, the role of J n cos(n0) may become relevant in the mixed state close to T c , for 
current I ~ I c and voltage < V < V c . 
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Intrinsic Josephson junctions are naturally formed in cuprates and correspond to pairs of 
Cu02 planes, separated by insulating layers. Such arrays have become attractive for quan- 
tum computation^^. Higher harmonics alter the profile of the tilted washboard potential 
and, thus, the sequence of tunnelling and dissipation processes determining the quantum 
device operation. Furthermore, arrays of intrinsic Josephson junctions biased in the resistive 
state where the Josephson current oscillates are being deployed as terahertz emitters^— . 
Upon decreasing the bias from the fully resistive state, the emission power increases as the 
Josephson frequency resonates with the cavity. With further voltage decreasing, some junc- 
tions may fall into the superconducting state, thus increasing voltage on the other junctions 
and, ultimately, switching off the radiation. Fine tuning and control of such intertwined 
oscillating-dissipative processes is crucial for the correct operation of the emitter. 

The critical temperature anomalies, observed when the ]?-doping of the CuC>2 planes is 
varied, are quantified. Such anomalies have been related to the emergence of a striped high- 
T c phase, with spatially modulated superconducting order, depending on the doping level 
p. An effective higher-order Josephson coupling varying as a cosine function of twice the 
difference of the superconducting phases on adjacent planes has been demonstrated. Several 
concomitant evidences of antiphase ordering in cuprates, besides the strong suppression of 
T c , have been reported^— . The T c suppression will be modeled by using the nonsinusoidal 
Josephson junctions with the ratio of the second to first harmonics depending on the doping 
level p. We remark that an array of junctions with simple sinusoidal current-phase, while 
correctly describes homogeneous low-T c superconductors characterized by uniform positive 
Josephson coupling, seems quite inadequate for the complex phenomenology of strongly 
correlated high-T c cuprates. 

II. NONSINUSOIDAL RSCJ MODEL 

A two-dimensional array of Josephson junctions is sketched in Fig.[T](a). The bias current 
lb is injected to the left electrode and collected from the right electrode. Circles represent 
superconducting grains connected by weak-links (crosses). According to the RSCJ model, 
the current Uj flowing through each junction is: 




(3) 
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where and i?y are the shunt capacitance and resistance between grains i and j, Is,ij{(j>ij) 
is the Josephson current, is the Langevin fluctuation source. The voltage drop across 

the junction is given by: 

13 1 3 2e dt ' 1 ' 

with (pij the phase difference of the order parameters in the grains i and j. In the usual 
RSCJ model, Is t ij{4>ij) is a simple sinusoid, whereas in the present work the nonsinusoidal 
form given by Eq. ([2]) is considered. Therefore, the current Iy flowing through each junction 
connecting the grains i and j writes as: 

J ij = Ci i~Jr + + ^ ^ n ' ij sin ( n ^i) + Jn,ij cos(n0y)] + SI L;ij . (5) 

Rij n >i 

Jjj is given by the sum of the following contributions: the charging current through the 
shunt capacitance CV,-, the Ohmic current through the shunt resistance R^, the n Josephson 
current sources I n>i j sin (ncpij) and J n>i j cos(n0y) and the Langevin current. 

It is worth noting that for n = 1, sin0jj is the familiar sinusoidal Josephson current 
I Ci ij sin0y, whereas Ji^j cos(n0y) with oc Vij/Rij corresponds to the voltage times a 
phase-dependent conductance term 1/ cos(0 i j)^ i ^ r — . Therefore, Eq. (E]) can be rewritten 
as: 

ft c% 1 ft c% . 

where Eq. (J4]) has been used. The second term on the right hand side of Eq. ([6]) is commonly 
called the interference current. 

The equivalent circuit of a junction obeying Eq. ([5]) is shown in Fig. Hlb). It corresponds 
to the parallel of a linear capacitor CV,-, a linear resistor R^, a parallel of n inductors L n ^j 
(related to the I n ^j sin(n0jj) terms) and a parallel of n memristors M n i j related to the 
Jn^j cos(n(pij) terms (we use the notation memristor after—). 

Eq. fl5]) can be written more compactly as: 

A? = C ij -JT + "TT + l WM Si n (™0ij + <fio,n,ij) + SI L ,ij , (7) 

dt K n > 1 

with: 



Ic,n,ij \j I>n,ij ^n,ij i (8) 
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and: 

<fio,n,ij = arctan ( -J^-J . (9) 

Conventional Josephson junctions are usually classified in terms of the Stewart- 
McCumber parameter f3 c = t rc /tj with r RC = RC and Tj = § /27rI c R , as overdamped 
(j3 c <C 1), general (j3 c ~ 1) and underdamped (j3 c ^> 1). For the nonsinusoidal junction 
described by Eq. (JTJ), the definition of the Stewart-McCumber parameter can be generalized 
as follows: 

# = ^ , (10) 

with 

$ 

Tj* = — . (11) 

27T X/n 1 c,n,ijRo 

Eq. (J7J) can be numerically solved for an arbitrary number n of harmonics. Nonetheless, 
we will restrict our discussion to the following two cases relevant for the applications: 

h,ij{<f>ij) = h,ij sin(0y) + J ltij coa(<f>ij) , (12) 

and 

h,ij{<t>ij) = h,ij sin(0ij) + I 2 ,ij Bin(2^ y ) . (13) 

The scheme of the current- voltage characteristics of an underdamped (/?* ^> 1) Josephson 
junction obtained by solving Eq. (JTJ) is shown in Fig. |2j In particular, Fig. |2£a) refers to 
the simple sinusoid, Fig. Mb) refers to Is,ij(4>ij) given by Eq. fTl2|) and Fig. Hfc) refers to 
Is,ij(4>ij) given by Eq. (fTBl . The intermediate states are characterized by voltage drops in 
the range < V^- < V Cji j and current = I c ,n,ij- Upon current (voltage) decrease starting 
from the normal state, the behavior is always resistive, implying that the system reaches the 
superconductive ground state without exploring the intermediate states. For overdamped 
junctions (/?* 1), the intermediate states are characterized by voltage drop and current 
respectively in the range < < 2V c>i j and I c , n ,ij < hj < Ic,n,ij[2Vc,ij}- Upon increasing and 
decreasing the external drive, the current- voltage behavior is the same, hence no hysteresis is 
observed. In the general case (0* ~ 1), the I — V curve is partly hysteretic. Upon increasing 
the external drive, the intermediate states are characterized by a voltage drop in the range 
< Vij < V C) ij and current equal to I c ,n,ij- As the external drive decreases, the backward 
current lies slightly below the forward current. It is worthy of remarks that the capacitive 
effect is reduced with the nonsinusoidal current phase relation in comparison to the simple 
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sinusoidal case. 

As a final remark, we note that since the simulations are addressed at modeling the zero- 
frequency (time-asymptotic) response of a macroscopic array, the Langevin term does not 
affect the results and thus in the simulations the term Sl^^j can be set to zero. Nonetheless, 
we stress that the term Sl^^j has profound conceptual implications related to the microscopic 
random dissipation/tunneling events and the onset of decoherence according to the Caldeira- 
Leggett picture. The term SlL,ij plays a major role in the evaluation of current noise power 
spectra shape and amplitude^. 



A. Resistive transition in granular superconductors 

The resistive transition is modeled by using a disordered network of weak-links with non- 
sinusoidal current-phase relation. The network is routinely solved by a system of Kirchhoff 
equations by using Eq. ([7]) and Eq. ([2]) in the temperature range just below T c . The network 
is biased by constant current 1^. The following conditions will be used for the simulation: 

1. The superconductive ground state of each weak link is characterized by current Iij < 
mm{I cni j} = I^min and = 0. The symmetry is not broken, thus J n ^ vanishes. The 
conductance of the weak-links in the superconductive state is taken G 3> e 2 /h 

i.e. G is much greater than the quantum conductance e 2 /h. This condition guarantees 
the existence of the superconductive ground state. 

2. The intermediate states correspond to the coexistence of superconducting and normal 
domains. According to the two-fluid model, unpaired electrons coexist with paired 
electrons in the region of temperature close to T c respectively, with densities: 

rto(T)=n (£) (14) 



T J 



l - [ T. 



T \ 4 " 

(15) 



where n is the total density of normal electrons. The fraction n§ of superelectrons is 
characterized by critical current I c , n ,ij — In,ij- Conversely, the fraction of normal 
electrons has a finite value of J n jj and, thus, from Eq. (j7J), is characterized by critical 
current I c ,n,ij = \Jln,ij + J%,ij- The condition < 1^ < yj + J 2 ^ holds in the 
intermediate state. The conductance of the weak-links in the intermediate states varies 
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between G and G a = 1/R , as a function of temperature, according to the relative 
fraction of super to normal electrons. 

3. The normal state is achieved when the voltage Vij across the junction exceeds V Cij . 
The conductance of the weak-links is G D = 1/R Q - The current flowing through each 



The superconductor-insulator transition is simulated by solving the system of Kirch- 
hoff equations at varying temperature. The critical currents I n ^j and J n ^ are assumed to 
vary on temperature according to the linearized equations I n ^j = I , n ,ij (1 — T/T c y and 
Jn,ij = Jo,n,ij (1 _ T/T c y , where I ,n,ij and J ,n,ij are the lowest temperature values of Inl- 
and J n> ij and the exponent 7 is about 2 for high-T c superconductors. Hence, the criti- 
cal current I c ,n,ij depends on temperature according to I c ,n,ij = h ,n,ij (1 — T/T c y , with 



Ic ,n,ij — y!o,n,ij + ^o,n,ij- ^ n or der to take into account the disorder of the array, I Uii j and 
Jn t ij are taken as random variables, distributed according to Gaussian functions with mean 
values 7 0>n and J 0)Tl and standard deviations A/ 0jn = AJ 0tTl . 

By effect of the temperature increase and consequent reduction of the critical current, 
the weak-link with the lowest value of the critical current I c>n ,ij — Ic,min switches to the 
intermediate state and, then, becomes resistive when Vij > V c . The resistive transition of the 
first weak-link has the effect to set the value of the voltage drop across the other weak-links 
in the same layer. The result is the formation of a layer of weak-links either in the resistive 
or in the intermediate state. As temperature further increases, the critical current I c , n ,ij 
further decreases. More and more weak-links gradually switch from the superconductive to 
the intermediate state and then to the resistive state. The term j njij - acts by increasing the 
critical current value of the weak-link in the intermediate state in the layers undergoing the 
transition. It is worthy to remark that the increase of critical current is relative to the fraction 
of normal electrons in the mixed states. The onset of J n> ij cos(n0jj) is indeed triggered by 
the elementary resistive transition of the weak-link with the lowest critical current, since it 
is related to the partly broken pair-symmetry of the weak-links in the intermediate state. It 
has no effect on the links in the superconductive state, neither on those in the fully resistive 
state. 

Fig.[3]shows the curves of the resistive transitions obtained with the current-phase relation 
Is,ij{4*ij) = sin(0jj) + l2,ij sin(20^) for a two-dimensional 30 x 30 network. The curves 



weak-link satisfies: 1^ > y In,ij + Jn,ij- 





S 



correspond to different values of the term li^j- The values of the critical currents are 
Ii s ij = 1mA and ranging from to 1mA. The standard deviation of the critical currents 
is A/ Coj „ = 0.3mA. The effect due to sin(2</>y) corresponds to a shift of the transition 
towards higher or lower temperature depending on amplitude. As opposed to J n ^j cos(n0^), 
the term 7 2 ,r/ sin(20^) acts on the Josephson coupling and thus its effect is higher at the 
beginning of the transition and decreases as the fraction ris of superelectrons decreases. 

Fig. H] shows the curves of the resistive transitions obtained with current-phase relation 
Is } ij(4>ij) = Iijj sin(0jj) + J\ t ij cos(0jj) for a two-dimensional 30 x 30 network. The curves 
correspond to different values of the term J\^. The values of the critical currents are 
= 1mA and ranging from to 1mA. The standard deviation of the critical currents 
is AJ on = 0.3mA. Initially, the weak-links are in the superconductive state, thus the 
network resistance is negligible. As temperature increases, the weak-link with the lowest 
critical current switches to the intermediate state and then to the resistive state with the 
consequent onset of the term J n ^j cos(n0jj) and redistribution of the currents. One can 
notice that the curves overlap at the beginning of the transition, whereas become more 
separated when T — > T c , implying that the effect of the term J n ^j cos(n0jj) is more relevant 
as the transition approaches its end. The amplification of the J n ^j cos(n0jj) effect, as the 
resistance increases, means that J n ^j acts as modulation of the resistance. The modulation 
effect due to J n ^j can be noted at the level of each elementary transition step. Fig. Hfb) and 
(c) show the zoom of the resistance steps corresponding respectively to the beginning and to 
the end of the transition curves in Fig. H(a). One can notice that the microscopic deviations 
from the staircase profile obtained with the simple sinusoidal current-phase relation increase 
as the global network resistance increases in agreement with the modulating action of the 
term J n>ij . 

B. Critical temperature anomaly in cuprates. 

In the previous section, the approach has been applied to a granular superconductor where 
the disorder of the material is taken into account by using a suitable probability distribution 
function of a relevant parameter. In particular, the probability distribution function is a 
Gaussian with the variance a accounting for the randomness of the critical currents I c over 
the array. In this section, the approach is implemented to model perfectly ordered single 
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crystals. In the absence of localized and extended defects, the relevant parameters of the 
arrays of Josephson junction are expected to be deterministic and thus a probability distri- 
bution function is not necessary, i.e. a = 0. 

Specifically, the proposed method will be applied to the intrinsically formed arrays of Joseph- 
son junction in single crystal of layered cuprates. The reported simulations are mainly ad- 
dressed to describe the occurrence of a 7r-phase shift within the CuC>2 planes in accordance 
with the anti-phase ordering model put forward by Berg et. al.— to account for many 
anomalies exhibited by cuprates. The existence of an anti-phase ordering has been exper- 
imentally confirmed The approach presented in this work is particularly suitable 
to simulate the array of intrinsic Josephson junctions with unconventional current voltage 
characteristics and the onset of 7r phase. Thus, the goal is a phenomenological descrip- 
tion of the model^ and the corresponding simulation of the experimental results presented 
j n 67- - 70_ j n p ar ticular, the predictions are compared with the experimental data concerning 
the T c suppression observed in cuprates at varying levels of doping. The doping p, i.e. the 
number of holes per copper atom in the CuC>2 planes, is a key quantity determining the 
main properties of high-T c superconductors, whose typical structure is shown in Fig. [5j A 
parabolic relationship between superconducting transition temperature T c and doping p has 
been envisaged: 

l--^- = 82.6(p-0.16) 2 . (16) 

-L c,max 

This relation is inaccurate for certain values of the doping and a very pronounced T c sup- 
pression (as shown in Fig. [6]) has been reported in many cuprates^— . Such a universal 
suppression of T c has been ascribed to the tendency of charge stripe formation, with spa- 
tially modulated superconducting order and phased. Evidence of stripe order in cuprates is 
provided by the enhancement of the anisotropy of resistivity with temperature. The charge 
dynamics is those of a superconductor in plane at high temperature but the behavior is that 
of a poor metal in the orthogonal direction. At low temperature, the effective Josephson 
coupling is always positive yielding the homogeneous low-T c superconducting phase, while 
the striped superconducting phase is found at relatively higher temperature. The scheme 
describing such phenomenon is given by alternating stripes of superconductor and insulator, 
forming an array of Josephson junctions. For a <i-wave superconductor at high temperature, 
with a strong crystal field coupling which locks the lobes of the pair-wave function along 
the crystallographic axis, the order parameter may change sign under rotation across the 
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planes. The overall result is an effective higher order Josephson coupling depending on the 
cosine of twice the difference of the superconducting phases on neighboring planes and the 
dominance of a negative sin (20) component in the current-phase relation^. 
Such an array of Josephson junctions, with spatially modulated Josephson coupling, exhibit- 
ing anomalous transport and thermodynamics, could not be accounted for by the resistively 
and capacitively shunted model with simple sinusoidal current-phase relation. Therefore, 
in the present work, a simulation based on arrays of nonsinusoidal Josephson junctions is 
put forward. The scope is the estimation of the critical current components yielding the 
parabolic dependence of T c and the suppression observed at doping values close to p — 1/8. 
The different doping level of the CuC>2 planes is taken into account by varying the critical 
current, which is related to the number of Cooper pairs in the superconductive phase and 
thus enhanced/suppressed by the hole doping. Specifically, the variation of the negative 
component 7 2 ,2j sin(20jj) dominating over the simple sinusoid is taken into account as the 
origin of the suppression of the critical temperature. As already stated, the simulations refer 
to a perfect crystal lattice instead of a granular superconductor, thus the array is perfectly 
ordered and a is negligible (a = 0). Hence, one can expect that the superconductivity 
is easily suppressed by current perpendicular to the superconducting layers while current 
flowing parallel to the layers would not destroy the superconducting state of the crystal. In 
the following, two sets of simulations are performed. 

First, the transition is simulated to obtain the ideal parabolic dependence given by Eq. (fT6l) . 
The differential equation ( Tl2|) is solved for nonsinusoidal junctions with and 
components. The I-V characteristics of the single junctions is obtained and then imple- 
mented to simulate the transition of the whole array as described in the previous sections. 
The transition curve allows one to deduce the critical temperature by using the relation 
In,ij — Io,n.ij (1 — T /T c ) 7 , where the exponent 7 is taken equal to 2. In the present work, the 
critical currents have been varied in the range 0.1 4- 10mA corresponding to doping level p 
varying between 0.05 -4 0.18 and critical temperature varying between -4 95K according to 
the data of Refs.— ~— . The critical temperatures and currents obtained from the simulation 
are shown in Fig. [7] (a) and (b) (blue squares). 

Then, the transition is simulated to obtain the suppression of T c with respect to the parabolic 
dependence. The suppression of T c is obtained by a decrease of as a function of the 
doping level for values ranging from p = 0.08 to p = 0.17. The critical temperatures and 
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currents obtained from the simulations are shown in Fig. [7J (a) and (b) (magenta circles). 
In the inset of Fig. [7(b) the ratio a of the second harmonics for the ideal parabola and the 
real curve with suppression, is plotted. The maximum temperature suppression corresponds 
to a value of the ratio close to 3.5. 



III. CONCLUSIONS 



The nonsinusoidal current-phase relation has been considered in the resistively shunted 
Josephson junction model for describing the superconductive transition. By solving a system 
of Kirchhoff equations for the array of nonsinusoidal Josephson junctions, it is found that 
additional cosine and sine terms modify the transition curves by changing resistance and 
Josephson coupling in the framework of the two-fluid model of superconductivity. Higher 
harmonics, besides the simple sinusoid I c sin 0, might arise in the vicinity of the transition 
because of the nonstationary conditions and the random distribution of potential drops and 
impurities in granular superconductors. The approach has been implemented for charac- 
terizing the critical temperature suppression observed in cuprates. In particular, our focus 
is on the anomalies experimentally observed in cuprates that need to go beyond the sim- 
ple sinusoidal picture arising from a constant positive Josephson coupling valid for low-T c 
superconductors. The specific example of the T c suppression at doping level p = 1/8 is 
described in terms of the ratio of the second to first sinusoidal components of the current- 
phase relation. The naturally formed networks of Josephson junctions, due to insulating 
layers sandwiched between Cu02 planes, are strongly affected by the presence of higher- 
order terms in the current phase relation. Further applications of the present approach can 
be envisaged to account for the complex phenomenology of high-T c materials forming arrays 
of nonsinusoidal Josephson junctions and its implications in novel quantum devices. 
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FIG. 1: (a) Two-dimensional Josephson junction array representing a granular superconductor. 
Circles represent superconducting grains. Crosses represent weak-links between grains. The bias 
current 1^ is injected to the left electrode and collected from the right electrode, (b) Equivalent 
circuit of the weak-link between the grains i and j. The linear resistor Rij, the linear capacitor 
Cij, the nonlinear inductor L n ^j and memristor M n ^j are connected in parallel. The current in- 
flows from grain i to grain j. Vij is the voltage drop across the weak-link. 
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FIG. 2: Josephson junction characteristics of a weak-link with current-phase relation (a) Is((/>) = 
I c sm(cp); (b) Is(<fi) = h sin(^>) + J\ cos(^) with I\ = 1mA and J\ = 0.5mA; (c) Is(c/>) = I\ sin(^) + 
I2 sin(2(/)) with I\ = 1mA and 1% = 0.5mA. The generalized Stewart-McCumber parameter is 
/3* = 45. 
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FIG. 3: Resistive transition of a two-dimensional network with current-phase relation of the form 
Is,ij{<t>ij) = sin(0jj) + l2,ij sin(2(^jj). The average value of the critical current is 1mA. 
The curves correspond to different average values of the critical current l2,ij, namely I^ij = 0mA, 
^2,ij = 0.5mA, l2 t ij = 0.75mA and i^i* = 1mA. The normal resistance R Q is 1Q equal for all the 
junctions. 
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FIG. 4: Resistive transition of a two-dimensional network with current-phase relation of the form 
Is,ij(4>ij) = h,ij sin(^>jj) + Jijj cos(4>ij). The average value of the critical current is 1mA. 
The curves correspond to different average values of the critical current Juj, namely J\m = 0mA, 
Ji.ij = 0.5mA, Ji ij = 0.75mA and J\^j = 1mA. The normal resistance R Q is lf2 equal for all the 
junctions. Panels (b) and (c) show the details of the beginning and the end of the transition. 
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FIG. 5: Arrays of intrinsic Josephson junctions are naturally formed in cuprates by the CuO"2 planes 
separated by layers of insulating atoms. The hole doping p of the CuO"2 planes affects transport and 
thermodynamic properties of cuprates. Several transport anomalies have been observed around 
p = 1/8 that cannot be explained in the framework of a conventional picture of the intrinsic 
Josephson junctions and have been ascribed to the antiphase ordering across the planes [52-58]. 
The modulation of the phase can be taken into account by using the proposed array of resistively 
and capacitively nonsinusoidal Josephson junctions. 
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FIG. 6: Critical temperature T c as a function of the hole doping p. The ideal parabolic relation 
is plotted as a reference (solid line). Circles are experimental data obtained on YBCO samples 
with varying doping level of the Cu02 planes (— ). The suppression of T c in the range of doping 
between 0.08 and 0.17 can be observed. 




22 



100 




0.00 0.05 0.10 0.15 0.20 



P 

FIG. 7: Critical temperature T c (a) and currents I c (b) as a function of the hole doping p. The 
data are obtained by simulating the network of Josephson junctions with current-phase relation 
given by _ + / 2 „ jS i„(2^. The average critical current of the array takes 

values in the range 1 mA 10mA. In order to yield the suppression of T c as a function of the 



